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We introduce a new family of networks, the Apollonian networks, that are simultaneously scale- 
free, small world, Euclidean, space-filling and matching graphs. These networks have a wide range 
of applications ranging from the description of force chains in polydisperse granular packings and 
geometry of fully fragmented porous media, to hierarchical road systems and area-covering electrical 
supply networks. Some of the properties of these networks, namely, the connectivity exponent, the 
clustering coefficient, and the shortest path are calculated and found to be particularly rich. The 
percolation, the electrical conduction and the Ising models on such networks are also studied and 
found to be quite peculiar. Consequences for applications are also discussed. 

PACS numbers: 89.75.Hc, 89.75.Da, 89.20.Hh, 64.60.Ak, 02.50.Cw, 05.50.+q 
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The study of networks has fruitfully inspired the un- 
derstanding of transport and information flow within sys- 
tems of many degrees of freedom during the last years 
[H Q, y, Ql . Many network types have been proposed and 
investigated serving to differentiate and elucidate effects 
ranging from clan-formation, epidemics spreading and lo- 
gistic planning to earthquakeprediction, neural activity 
and immunological defenses [l|, Q, y, |j{ . 

Among a multitude of models, much attention has been 
dedicated recently to the study of scale-free networks, 
namely, networks that display a power-law degree distri- 
bution, p{k) ex fc"^, where k is connectivity (degree) [l|. 
The widespread appearance of networks with fat-tailed 
degree distribution in Nature has been mainly justified 
by the fact that their evolution is self-organized in a way 
that hub elements are naturally generated and represent 
the dominant constituents of the network |5| . 

There are, however, several conceptual and applied 
aspects of scale-free network models that are still un- 
der investigation. For instance, the fact that networks 
with exponent 7 > 2 can be embedded in regular Eu- 
clidean lattices has motivated Rozenfeld et at. 6] to de- 
velop a method for generating random scale- free networks 
with well defined geographical properties. The ultrasmall 
characteristic found for the diameter of random scale- free 
networks -J] is also an intriguing aspect of these struc- 
tures that is expected to explain anomalies in their dif- 
fusional and transport phenomena behavior. Finally, the 
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stability to crashes due to removal of sites in such net- 
works has also been the focus of intense research activity 
;.8j. These studies revealed that random scale- free net- 
works remain connected even if nearly all of their nodes 
break down, as long as the degree exponent 7 < 3. As 
a consequence of this striking property, the Internet and 
other well known networks displaying scale-free behavior 
must be resilient to random breakdown. 

An important issue that has been much less explored 
in network science is the design problem of a suitable 
complex network for a given system where connectivity 
plays an essential role. For this purpose, deterministic 
networks represent the ideal strategy, since they are con- 
trollable and can therefore be subjected to the particular 
constraints of the problem, for example, those related to 
space, time, economical and other types of limitations. 
Using deterministic models, it should be possible to out- 
line some specific guidelines for the design of real net- 
works with desired properties (e.g., high clustering, com- 
pactness, etc.). Some examples of deterministic mod- 
els that are scale-free have been recently proposed and 
successfully used to describe random growing networks 
0, [i3 , but they have eluded the possibility of being em- 
bedded in Euclidean space. In this letter we introduce a 
class of networks, named Apollonian networks, that can 
be either deterministic or random, are scale-free, display 
small world effect, can be embedded in an Euclidean lat- 
tice, and show space-filling as well as matching graph 
properties. 

To define our network in its simplest deterministic ver- 
sion let us start with the problem of space-filling packing 
of spheres according to the ancient Greek mathematician 
ApoUonius of Perga [llj. In its classical solution, three 




FIG. 1: (a) Classical Apollonian packing, (b) Apollonian 
network; first, second and third generation are symbols Q, ■ 
and • respectively. 
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FIG. 2: Shortest path I as a function of the number of sites 
A^. The solid line is a guide to the eye with slope 3/4. The 
inset shows the semi-log plot of the clustering coefficient C 
as a function of N . The dashed line is a guide to the eye for 
C = 0.828, observed at large values of N . 



circles touch each other and the hole between them is 
fiUed by the circle that touches all three, forming again 
three but much smaller holes that are then filled again 
each in the same way as shown in Fig. QJa). The circle 
size distribution follows a power-law with exponent of 
about 1.3 |lj. In fact there exist many other topologies 
with different fractal dimensions [iJI . This procedure can 
also be generalized to higher dimensions T3] . 

Apollonian tilings can describe dense granular pack- 
ings and have also been used as toy models for turbulence 
and fragmentation. Connecting the centers of touching 
spheres by lines one obtains a network which in the clas- 
sical case of Fig. ^a) gives a triangulation that physically 
corresponds to the force network of the packing. We will 
call the resulting network shown in Fig. ^b) an "Apollo- 
nian network" . Besides resembling the graphs introduced 
by Dodds |lj| for the case of random packings, this net- 
work has also been used in the context of porous media 

First we show that Apollonian networks are scale-free. 



FIG. 3: Conductance G as function of the number A'^ of sites 
when a potential is applied either between the center and all 
the sites of the last generation (upper curve) or between two 
corners of the network (lower curve). The inset shows the 
distribution of local currents for the first case. 



From Fig.^b) one can see that at each generation n the 
number of sites N increases by a factor three and the 
coordination of each site by a factor two. More precisely, 
at generation n (n = 0, 1, 2,...) there are m{k,n) = 3", 
gra-i^ 3"~^..., 3^, 3, 1, and 3 vertices with degree k = 3, 
3 X 2, 3 X 2^,..., 3 X 2"-i, 3 x 2", and 2"+\ respectively, 
where the last number of vertices and degree correspond 
to the three corners, Pi, P^, and i-3. Due to the discrete- 
ness of this degree spectrum, it is convenient to work with 
the cumulative distribution P{k) — X]fc'>fe '^(^'j'^)/-'^"! 
where Nn = 3-1- (3"+^ — l)/2 is the number of sites at 
generation n. In order to compute the analog to our 
problem of the degree exponent 7, as commonly defined 
for continuum distributions i1DI|, it is easy to show that 
P{k) ock^-'', with 7 = l + ln3/ln2w 2.585. 

Next, we show that Apollonian networks display small- 
world effect ;3i, li ilQj . Accordingly, this means that the 
average length of the shortest path I between two vertices 
as defined in ref. u7| should grow slower than any positive 
power of the system size N. For the network in Fig.^b) 
we find that I oc (InA^)^, with /3 « 3/4 as can be seen 
in Fig. 13 This is a novel intermediate behavior between 
small {I oc InA'^) and ultrasm,all (I oc InlniV) networks 
i. 

The second property characterizing a small world net- 
work is the clustering coefficient C as defined in ref. [i3 ■ 
We find it to be 0.828 in the limit of large iV, as shown in 
the inset of Fig. [3 This is a large value, comparable for 
instance with the cluster coefficient for the collaboration 
network of movie actors C = 0.79 [l3. Consequently, 
since L grows logarithmically and C is close to unity, 
the Apollonian network indeed describes a small world 
scenario p^ . 

As opposed to hierarchical lattices [lO| used for in- 
stance in the Migdal-Kadanoff approximation, Apollo- 
nian networks are embedded in Euclidean space and can 



therefore describe geographical situations. In the case of 
Fig. nja) the sites form a multifractal point set in the 
asymptotic limit of infinitely many iterations [T^ . In the 
Apollonian networks for n — > oo, the bonds do not only 
completely cover the space (like a Peano curve), but also 
never cross each other. In this situation we have what is 
called a "matching graph" which is not the case for any 
other known scale-free network 1] . 

Now let us explore the properties of some well-known 
simple models defined on the lattice shown in Fig. ^b). 
Thinking of an electrical supply system where the central 
site is the source (power station) , the last generation sites 
are sinks (consumers), and all the connections have the 
same conductance, one can calculate the equivalent con- 
ductance G of the network. Going from one generation 
to the next, the number of parallel system outputs in- 
creases by a factor three while at each node the number 
of outputs doubles. Due to the self-similarity, one has 
therefore doubled the conductance while increasing the 
number of sites by three. As a result, G follows a power- 
law, G (X N^, with z — 2/3. As shown in Fig. O this 
exponent can be confirmed through computational sim- 
ulation of the resistor network model. From Kirchhoff 's 
law, the resulting system of linear algebraic equations 
are numerically solved using a sparse matrix inversion 
technique. This solution also allows us to compute the 
distribution of local currents P{I), as shown in the inset 
of Fig. O Interestingly it follows a power-law for large 
values of / with an exponent —1.24, as a consequence of 
the hierarchical choice of the output sites. 

One can alternatively apply a potential drop between 
two corners of the system, i.e. Pi and P2 in Fig. ^h) . 
The resulting conductance numerically calculated as a 
function of the size of the network gives a power-law too, 
but with exponent 0.45 as shown in Fig. |3| 

We now proceed with the study of bond percolation si- 
multaneously applied among the three points Pi , P2 and 
P3 of Fig. n^b) removing, however, the three direct con- 
nections between these points. Similarly to ref. Q, we 
find that the percolation threshold Pc goes to zero as the 

number of sites increases, according to a power-law, as 

_ j_ 

seen in Fig. 0| Assuming pc oc Vn " , the slope yields 
an exponent v fa S. 

Spin models can be defined on the Apollonian lattice 
in a variety of ways. We consider here that an Ising 
spin Si is placed on each site of the lattice. Between two 
sites we have a pair interaction constant J„, where n indi- 
cates that the link between the interacting spins is placed 
at the n-th generation of the construction of the lattice. 
Moreover, the interactions are considered to be ferromag- 
netic Jn > 0, decaying with n as J„ oc n~", and insensi- 
tive to the actual distance between the spins. With this 
model definition, it is straightforward to use a transfer 
matrix formalism to write down temperature dependent 
recurrence relations for the free energy /„ and correla- 
tion length ^n of the model at two subsequent genera- 

tioIlSj/n+i = fn+l{fn,£,n;T) and Cn+l = ^„+l (/„, ^n; P) 

|19l I20I . Exact values for these and any other functions 
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FIG. 4: Double logarithmic plot of the percolation threshold 
Pc as function of the inverse system size A''. 



in the thcrmodynamical limit n — > 00 can be obtained 
through the numerical iteration of fn+i and ^n-i-i and for 
their derivatives. The iteration process n is stopped when 
the free energy per spin is the same in two subsequent 
generations within the first 16 digits. 

The results obtained here are quite different from those 
for other self-similar lattices. The free energy, entropy, 
specific heat and magnetization (respectively f,s,c,m) 
are found to be completely smooth for all values of T. A 
classical Schottky maximum dominates the behavior of c 
while s increases monotonically to /jBln2. The presence 
of long-range order is observed as m = 1 at low values 
of T. On the other hand, when T is large, m ^ as 
exp(— T) for a = 0, with no noticeable presence of a crit- 
ical temperature. For a > 0, we find a stronger decay, as 
m ^ exp(— T"^). The results for the correlation length are 
more subtle. ^ diverges at low temperatures, expressing 
the evidence of long-range order. As T increases beyond 
a given Tc{n) it converges to a well defined value. How- 
ever, if the iteration procedure is pursued to a value of 
n larger than the one required for the convergence of /„, 
we observe that ^ keeps on diverging within a larger T 
interval. This behavior therefore suggests the existence 
of an n dependent effective critical temperature Tc{n), as 
has been found by other authors for spin models on an- 
other scale-free lattice |2lj. In Fig. we show how Tc{n) 
depends on n for different values of a. The power-law 
Tc{n) ^ n'^^°'\ with r going continuously from t(0) = 1 
to t(cxd) — 1/2 indicates a novel kind of critical behavior. 

It is important to mention that a more general version 
of the network shown in Fig. ^h) is to choose at genera- 
tion n in each triangle a new site randomly and connect 
it to the three adjacent sites. In this sense the network 
can describe the porous fine structure inside a block of 
heterogeneous soil surrounded by larger pores or as the 
network of roads or supply lines going to smaller urban 
units surrounded by larger ones. There exists in fact a 
very large class of Apollonian networks obtained not only 
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FIG. 5: Double logarithmic plot of the critical temperature 
as function of the generation number n for different exponents 
a. 



by changing the choice of the position of the sites inside 
the areas or changing the initial graph but also having a 
different topology as classified in ref . J22| . All these net- 
works have the properties presented here, however, with 
different exponents if the topology is different. 

Apollonian networks have the typical properties of the 
complex graphs in nature and society: large clustering co- 
efficients, slower than logarithmic increase of the shortest 
path with number of sites, and a power-law degree distri- 
bution. In the context of critical phenomena, we did not 
find a finite critical point neither for percolation nor for 
the Ising model on Apollonian networks. The exponent 
1/ — 3 controlling how the percolation threshold goes to 
zero with N is unusually high while Tc{n) of the Ising 
model diverges with a power-law in the number n of gen- 



erations. 

Considering the description of the force network in a 
polydisperse packing, the small world nature of the lat- 
tice implies that between any two grains there are only 
a few contacts at best. If the medium is preconstrained, 
the electric conductance and elastic response have the 
same divergence |23|, and we find that the system be- 
comes more rigid with the inclusion of smaller particles 
following a power-law (with exponent 0.45 in our topol- 
ogy) in the number of particles. Giving to each grain a 
remnant magnetisation as discussed in ref. J24| . our re- 
sult for the Ising model implies that, in the limit of many 
particles, one always find a ferromagnetic response. 

Interpreting our network as a porous medium, the fact 
that Pc = is compatible with Archie's law, namely, that 
there exists no finite critical porosity and the medium 
always percolates. If the network represents roads or 
electric supply lines between cities implies that they have 
small world character and that the distribution of car flux 
or electrical currents follows a power-law when going from 
larger to smaller sites. The results for the Ising model 
imply that opinions or cooperations can stabilize at any 
degree of interaction, and the percolation threshold being 
zero means that in the thermodynamic limit the global 
connectivity cannot be disrupted, no matter how many 
connections are broken. 

We have seen that the Apollonian networks introduced 
in this paper can have several applications yielding rea- 
sonable conclusions while on the other hand their criti- 
cal properties are rather special. They should be stud- 
ied in more detail including synchronization, spreading, 
path-optimization, search algorithms and more specific 
features of the real systems. 
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